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Thon rings from amorphous ice and implications of beam-induced 
Brownian motion in single particle electron cryo-microscopy 
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Abstract 

We have recorded dose-fractionated electron cryo-microscope images of thin films of pure flash-frozen amor¬ 
phous ice and pre-irradiated amorphous carbon on a Falcon II direct electron detector using 300 keV elec¬ 
trons. We observe Thon rings [1] in both the power spectrum of the summed frames and the sum of power 
spectra from the individual frames. The Thon rings from amorphous carbon images are always more visible 
in the power spectrum of the summed frames whereas those of amorphous ice are more visible in the sum of 
power spectra from the individual frames. This difference indicates that while pre-irradiated carbon behaves 
like a solid during the exposure, amorphous ice behaves like a fluid with the individual water molecules un¬ 
dergoing beam-induced motion. Using the measured variation in the power spectra amplitude with number 
of electrons per image we deduce that water molecules are randomly displaced by mean squared distance of 
~ 1.1 A 2 for every incident 300 keV e _ /A 2 . The induced motion leads to an optimal exposure with 300 keV 
electrons of 4.0 e _ /A 2 per image with which to see Thon rings centred around the strong 3.7A scattering 
peak from amorphous ice. The beam-induced movement of the water molecules generates pseudo-Brownian 
motion of embedded macromolecules. The resulting blurring of single particle images contributes an addi¬ 
tional term, on top of that from radiation damage, to the minimum achievable B-factor for macromolecular 
structure determination. 
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1. Introduction 


2. Theory 


During the 1980s, Dubochet and his colleagues 
[2U3] developed a method for carrying out electron 
cryo-microscopy (cryoEM) of biological structures 
embedded in thin films of amorphous ice. Their 
early work involved careful comparison of the con¬ 
ditions needed to obtain thin films of amorphous 
rather than hexagonal or cubic crystalline ice [2]. 
The method they developed consisted of rapidly 
freezing a thin film of water, or buffer solution, 
containing the biological structure of interest by 
plunging it into liquid ethane at a temperature just 
above that of liquid nitrogen. Images of the result¬ 
ing specimen recorded using a suitable microscope 
could then be analysed to determine the structure. 
The method has since grown immensely in its power 
and popularity as improvements in technology have 
transformed its capability [3]. 

There has long been a debate within the electron 
cryo-microscopy community as to whether Thon 
rings [T] can or should be observed in power spec¬ 
tra of images of pure plunge-frozen amorphous ice. 
Thon rings arise naturally in bright-held phase con¬ 
trast images of any amorphous material and their 
origin is described in most standard texts 13 El- 
The spacing and shape of the rings depend on the 
electron-optical parameters that describe the im¬ 
age, principally the amount of defocus and astig¬ 
matism. Thon rings from amorphous carbon are 
routinely used to adjust astigmatism and set the 
defocus of a microscope. Amorphous ice has ap¬ 
proximately one half the density of amorphous car¬ 
bon but as oxygen atoms scatter more strongly than 
carbon atoms [5] one might naively expect to see 
Thon rings of similar strength from amorphous ice 
and amorphous carbon. 

In this paper we show that with the recently in¬ 
troduced CMOS direct electron detectors it is pos¬ 
sible to see unambiguously Thon rings in images of 
high purity amorphous ice. The higher detective 
quantum efficiency, DQE, of the new detectors al¬ 
lows weaker signals to be seen and the ability to col¬ 
lect images continuously in a dose fractionated (or 
movie) mode allows optimal exposure conditions to 
be chosen and new modes of image processing to be 
used, long after the sample has been removed from 
the microscope. 


2.1. Definitions 

In this work the pixel value in an image at pixel 
(r, s ) is denoted f r , a , and scaled to units of e _ /pixel. 
The discrete Fourier transform, /(it), at spatial fre¬ 
quency it is calculated using 

f(u) = ^2 frs exp(—i27ru-av s ). (1) 

r,s 


and the power spectra, S(u), calculated using the 
convention 


S{u) = -r^~\f(u)\ 2 = 1 


N u 


Nt, 


T, frs exp(—i2irx rs .u) 
r,s 

( 2 ) 


2.2. Noise whitened power spectra 

The modulation transfer function, MTF, and 
detective quantum efficiency, DQE, of the Falcon II 
detector are given in jTj. The MTF of the Falcon II 
drops rapidly with increasing spatial frequency cor¬ 
responding to the Falcon II having a broad point 
spread function, PSF. For the combination of pixel 
size and sensitive layer thickness used in the Fal¬ 
con II the PSF is chiefly determined by charge car¬ 
rier diffusion in the sensitive layer rather than the 
scattering of incident electrons [7,. As a result the 
DQE of the Falcon II does not fall with the MTF. 
The strong spatial frequency dependence of the MTF 
is however reflected in the power spectra from the 
Falcon II and can make it difficult to see weak sig¬ 
nals. In this case it is advantageous to use the noise 
whitened power spectrum, VF(it), given by 

W{u) = S(u)/N{u) (3) 

in which jV(tt) is the normalise noise power spec¬ 
trum [7]. For uniformly illuminated images, noise 
whitening removes the spatial frequency dependence 
of S(u) by boosting the signal at higher spatial fre¬ 
quencies. In particular it gives a flat background 
on which small signals can more easily be seen. 

In [7] it was shown that for the Falcon II, Af(u) 
could be estimated directly from the measured mod¬ 
ulation transfer function, MTF. Using J\f(u) esti¬ 
mated in this way results in only a few percent 
variation in W (it) over the range from near zero 
to the Nyquist frequency in images with no sam¬ 
ple present. In the present work it was found that 
the residual variation in W (it) with it could be fur¬ 
ther reduced by using a fit to the measured S(u) 
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with model for J\f(u) consisting of only a radially 
symmetric function and its low order aliased terms. 
Note that while A f(u) is a two dimensional function 
almost all of the non-circularly symmetric compo¬ 
nents arise simply from the contributions of aliased 
terms. Details of the fitting proceudure are given 
in the supplementary data. 


2-4 ■ Power spectra of Dose fractionated series 
If a dose fractionated exposure consists of M 
frames with an average of d uncorrelated electrons 
per frame, from Eq. |7| the expected value of Wm,™ 
defined as the sum of the M/m noise whitened 
power spectra from the images obtained by sum¬ 
ming to consecutive frames is 


2.3. Probability distribution for values in the noise 
whitened power spectra 

For uniformly illuminated images the probabil¬ 
ity distribution for values of W(u) can be described 
using a single parameter. With no sample, the real 
and imaginary components of /(it) in Eq. ([I]) are 
independent Gaussian random variables [8j. As a 
result the squares of the real and imaginary compo¬ 
nents have \ 2 distributions and their sum, S(u), an 
exponential distribution. Noise whitening sets the 
same average value, and hence exponential distribu¬ 
tion for values of W (it), at all spatial frequencies, it. 
For any given spatial frequency, it, the probability 
that W{u) has value y is therefore 

Prob[y = W{u)\ = exp(— y/T)/T (4) 

in which T is the mean value of W(u). 

The DQE of the Falcon II detector is principally 
limited by the intrinsic variability of its response to 
individual electrons. If the average gain and vari¬ 
ance in the gain are g and cr^, respectively, then 
the value of noise power spectrum in the limit of 
zero spatial frequency with n incident independent 
electrons per pixel is [9] 


M/m 

W M , m = £ (md/DQE) = ^. (8) 

i ^ 

in which D to t = Md is the total number of elec¬ 
trons in the exposure. Similarly the expected noise, 
Nm ,m 1 in Wm,™ is 


N M.m — 


\ 


M/m 

Y (TOd/DQE) 2 


Aot frn 
DQE V M' 


(9) 


While the value of Wm,™ is independent of to, the 
noise in sum, Nm,™, grows as the square root of to. 
In practice residual correlation in pixel values be¬ 
tween frames, such as from an offset drift, will result 
in an to dependence in Wm,u- The actual Falcon 
II detector used here had a small correlated shift 
in the row-reset between successive frames that re¬ 
sults in a 1% drop in WM, m {u ) between to = 1 and 
to = 2. However the largest systematic correlation 
comes from the applied per pixel gain correction of 
the images which results in an increase in Wm,™ (u) 
with increasing to but with careful gain calibration 
keeps the effect is less than 5% in going from m = 1 
and to = 120. 


S(0) = g 2 n +cr 2 n (5) 

and 

DQE(O) = g 2 / (g 2 + a 2 ). (6) 

Using Eqns. (|3| to ([b]) gives 

W( 0) = W(u) = n/DQE(0) = T (7) 

since the average value of W (tt) is the same for all 
it, by definition Af{0) = 1 and expressing the de¬ 
tector output in units of incident electrons sets the 
gain to unity. The probability distribution for val¬ 
ues of W(u) depends solely on T = n/DQE(0) via 
Eq. Q and has both mean and standard deviation 
of E. In particular, for the Falcon II DQE(O) ~ 0.5 
[7| and so T ~ 2 n. 


2.5. Effect of random motion of water molecules on 
power spectra 

At any given instant there is typically only a sin¬ 
gle high energy electron interacting with the sam¬ 
ple during a typical cryoEM exposure. An incident 
electron passing through a thin layer of amorphous 
ice will see a particular conformation of atoms but 
the radiation damage resulting from the interaction 
of the high energy electron with the amorphous 
ice will perturb the conformation of atoms in the 
ice [TO]. Subsequent high energy electrons incident 
on the ice will see, and perturb, the atomic con¬ 
figurations resulting in a Markov-like process for 
the evolution of the atomic configuration. In real¬ 
ity possible atomic configurations in amorphous ice 
are strongly constrained by the ice rules of Bernal 
and Fowler (TTi. The simplest approximation for 


3 






the transition probability of atom positions in go¬ 
ing from one conformation to the next is to use a 
Gaussian distribution depending solely on the num¬ 
ber of high energy electrons incident per unit area 
between the configurations . In particular if after 
n electrons per unit area the position of the i-th 
water molecule is Xi{n), after a further d electrons 
per unit area the probability distribution for its po¬ 
sition has a Gaussian distribution 


d = D/M we have S(u) oc D 2 . On the other hand 
if there is no correlation between frames then a u —> 
oo and S(u) —» Md 2 F 0 (u) 2 . 

For an exposure with a total of D electrons but 
where the molecules are continuously moving we 
can let M —► oo and the summation in Eq. (Ill 
goes over to an integral giving 


S(u) = 2F 0 (u) 2 [a u D + exp(-a u D) - l]/a 2 (15) 


P{xi(n + d) : ®j(n)) = 


(27T(i<7g) 3 / 2 


exp s - 


| Xj(n + d) -Xj(n )| s 

2d<7Q 


( 10 ) 


in which ctq is the mean squared displacement in 
a given direction in response to a single incident 
electron per unit area. 

Assuming an image consists of M frames and 
has a total exposure of D electrons per unit area. 
The power spectrum, S(u ), at spatial frequency u 
in terms of the Fourier components, fi(u), of the 
i-th frame is 

M M 

S(u) OC <|^/j(M)| 2 > = <fi( U )fj( U )> 

i—1 i^j — 1 

( 11} 

where < ... > denotes an ensemble average. Since 
the average number of incident electrons per unit 


area in each frame is d = D/M , using Eq. (10) 


the |* — j\d incident electrons between i-th and j-th 
frames will randomly displace the water molecules 
so that 

<f*( u )fj( u )> « d 2 F 0 (u) 2 exp(-a u \i - j\d) (12) 

in which a u = 2n 2 OgU 2 and cr 2 is the induced mean 
squared motion of a water molecules per incident 
electron though a unit area while Fg(u ) 2 is defined 
as 


F 0 (u) 2 = 


1 


M 


Md 2 


fi(u) 




(13) 


Substituting Eq. (12) into Eq. (Ill and summing 
over i and j gives 

S(u) = d 2 F 0 (u) 2 x 

o p -a u d 

{(1-e-““ d )M + . 


(1 - e -““ d ) 2 


— 1} + M ). 

(14) 


which in the limits of small and large a u goes to 
D 2 F 0 (u) 2 and 2DF 0 (u) 2 /a u , respectively. 

2.6. Estimation of beam induced movement from 
dose fractionated images 
The amount of movement resulting from an in¬ 
cident electron can be estimated from the variation 
in Thon ring signal as a function of the number of 
incident electrons in an image. When the induced 
motion is smaller than a given spatial frequency 
the Thon ring modulation at that spatial frequency 
is expected to grow quadratically with the number 
of incident electrons in an exposure but as the in¬ 
duced motion becomes comparable with the spatial 
frequency the amplitude should grow more slowly. 


From Eq. (15) the power spectrum, S(u), at a spa¬ 
tial frequency u, is expected to vary with total num¬ 
ber of electrons per unit area, D, as 


S(u) oc 2 [a u D + exp (—a u D) — l]/a 2 


(16) 


. (Jq is the mean-squared 
movement in a given direction per incident elec¬ 
tron passing through a unit area. The actual Thon 
ring modulation does not appear explicitly in this 
equation but the Thon ring modulation enables the 
signal from the ice to be identified over the back¬ 
ground noise in the measured power spectra. The 
measured spatial frequency dependence of S(u) will 
also include that of the detector but for the Fal¬ 
con II detector this can be removed by using the 
noise whitened power spectra, W{u) as defined in 
Eq. (§. 

From a single dose fractionated exposure it is 
possible to estimate Ug through the dose depen¬ 
dence present in Eq. (161. If an the exposure con¬ 


sists of M frames each with an average of d elections 
per unit area, the result of summing the original 
frames in blocks of m is to produce M/m frames 
with on average md electrons per unit area. Denot¬ 
ing the sum of the resulting M/m noise whitened 


In cases where there is very little induced movement 
a u —► 0 so that S(u) —> M 2 d 2 F 0 (u) 2 and since 
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power spectra by WM,m(u) and using Eq. (16) gives 
WM,m{u ) oc -—— (a u md + exp(— a u md) — 1 )/a u md 


(17) 

in which D tot = Md. The behaviour of Eq. © 
with 77i depends only onz = a u md and is described 
by the functional form 

g(z) = (z + e~ z -l)/z. (18) 

This is a monotonically increasing function that ini¬ 
tially varies linearly with 2 but plateaus towards 1 
for large 2 . By fitting the non-linear behaviour of 
WM, m ( u ) as a function of m it is possible to esti¬ 
mate a ui and hence obtain ctq. 


3. Experimental 

Quantifoil Rl.2/1.3 Cu/Rh grids were used to 
make a thin film sample of double distilled water 
(18 MSlcm conductivity). The grids were glow- 
discharged in air for 30-40 seconds and 3 /d of water 
applied in an environmental plunge-freeze appara¬ 
tus [15] ■ Grids were blotted for 4-6 seconds and 
rapidly frozen in liquid ethane. The grids were 
transferred to Krios cartridges and imaged using 
a Falcon II direct electron detector on a FEI Titan 
Krios operated at 300 keV. A nominal magnifica¬ 
tion of 75,000x corresponding a calibrated value of 
134,600x was used which results in a 1.04 A sam¬ 
pling with the 14 fj,m pixels of the Falcon II de¬ 
tector. Imaging was performed using nano-probe 
mode with parallel illumination using a 70 gm C2 
aperture and no objective aperture. A beam slightly 
larger and centred on the Quantifoil hole was used. 
Astigmatism and beam tilt correction were performed 
at the imaging magnification. Images were recorded 
at 18 frames/sec for either 8 or 17 seconds with ei¬ 
ther 2.33 e _ /frame, or 0.85 e _ /frame, respectively. 
The electron exposure per frame was determined 
from the screen current and average signal from the 
Falcon detector, with these being calibrated against 
a picoammeter as described previously [7]. 

The amorphous carbon control specimen was 
made in an Edwards 301 vacuum coating unit by 
evaporating carbon from an arc onto mica. The 
thickness of the carbon film was estimated to be 
150 A from the optical density of the carbon de¬ 
posited on adjacent piece of filter paper. The car¬ 
bon film was floated off the mica using a water bath 


and placed onto a Cu/Rh grid. The continuous car¬ 
bon grid was examined using the Atlas component 
of the FEI EPU software. From the Atlas an area 
of the carbon film was selected in which the carbon 
film had no holes or wrinkles both in the grid square 
and surround grid squares. After pre-irradiation of 
the carbon, images were recorded for 1.5 seconds 
with the same magnification and number of inci¬ 
dent electrons per frames as used for the ice sample. 
The stage was moved to adjacent areas in a raster 
(using 0.4 /im steps in X followed by 0.4 gm steps 
in Y). A total of 25 images from different but ad¬ 
jacent areas of carbon film were collected with the 
defocus for all images set as close as possible to a 
nominal 5000 A by manually adjusting the focus at 
each position. Later evaluation using CTFFIND3 
m showed that the defocus was actually 5990 A, 
with 200 A of astigmatism and ±90 A defocus vari¬ 
ation over the 25 images. 

Images were recorded using the full 4kx4k out¬ 
put of the Falcon II detector with the individual 
frames of all the images being captured using an 
in-house data capture setup. In processing the ex¬ 
posures the individual frames were not aligned com¬ 
putationally. 


4. Results 

Fig.g shows two types of noise whitened power 
spectra, W(u), J7| from a dose fractionated expo¬ 
sure of pure amorphous ice consisting of 141 frames 
with on average 2.33 e _ /pixel per frame. The spec¬ 
trum from the sum of all the 141 frames is shown 
in Fig. [lja) and with such a relatively high dose 
(300 e~/A 2 ) it is possible to see faint Thon rings. 
Fig. Qb) shows the sum of the 141 noise whitened 
power spectra of the individual frames. Exactly the 
same images were used but in the sum of the power 
spectra of the individual images the Thon rings that 
are barely visible in Fig. 0a) are clearly visible and 
can be seen to extend out beyond 3.4 A resolution. 
The -\/141 reduction in the noise from averaging the 
141 individual spectra enables the strength of the 
Thon rings in individual frames to be seen clearly. 
The strength of this signal indicates that the Thon 
rings arise from the intrinsic bulk water and not 
from impurities absorbed from the adjacent carbon 
film or contamination occurring during the blotting 
and transfer steps. 

The behaviour of the circularly averaged power 
spectra and circularly averaged noise whitened power 
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Figure 1: Thon rings from a dose fractionated exposure (im¬ 
age #230804) of amorphous ice obtained using either (a) 
the NWPS of the sum of the frames, or (b) the sum of the 
NWPS from each frame. The exposure consisted 141 frames 
at a dose of 2.33 e — /pixel using 300 keV electrons, a defocus 
of 7700 A and a pixel sampling of 1.04 A. The edge of the 
transform is at 1/2.08 A -1 and the strong ring pattern in 
(b) is centred around 1/3.7 A -1 . In both (a) and (b) the 
the images were first scaled so that the RMS noise level was 
1 and then the lower and upper grey limits set at -0.7 and 
1.3 about the mean values. 


spectra are compared in Sec. |SVl| of the supplemen¬ 
tary material. The flat, featureless backgrounds in 
Fig-0 a ) and Fig. 0 b) illustrate the success of the 
noise whitening procedure. The measured average 
value and noise in the noise whitened power spec¬ 
tra of Fig . |l| are within 5% of the values predicted 
by Eqs . (0) and <§ using n = 2.33 and DQE(0)= 
0.5. The probability distribution based on all the 
values in the 141 noise whitened power spectra used 
in Fig. 0b) is given in Sec. SVII of the supplemen¬ 
tary material. The distribution is well described by 
Eq. Q with a measured value of T = 4.91 that is 
within 5% of an estimated based on Eq. 0 again 
using n = 2.33 and DQE(0)= 0.5. 

To illustrate the origins of the difference be¬ 
tween Fig. 0a) and Fig. 0b) we carried out two 
control experiments using amorphous carbon film. 
The amorphous carbon film was first pre-irradiated 
as even atoms in films of carbon prepared by evap¬ 
oration from a carbon arc in vacuo move when ini¬ 
tially irradiated. The amount of movement de¬ 



Figure 2: Comparison from a dose fractionated exposure of 
an area of pre-irradiated carbon of (a) the NWPS obtained 
from the sum of all the frames and (b) the sum of the NWPS 
of individual frames. The exposure consisted of 25 frames 
with on average 2.33 e — /pixel per frame recorded using 300 
keV electrons, a 5990 A defocus and a 1.04 A pixel sampling. 
The grey scales of the images were set as in Fig.J^ 


creases with exposure but ~100 e“/A 2 is sufficient 
to effectively stabilise a film, i.e., the observed Thon 
ring pattern showed no drift and was stable. As 300 
keV electrons can cause displacement damage [T3], 
there was some residual movement but relative to 
the initial movement this was negligible. 

In the first control experiment, the initial 25 
frames from a 1.5 second exposure of an area of 
pre-irradiated carbon were used. The magnification 
and number of incident electrons per frame were 
the same as in Fig. 0 The corresponding noise 
whitened power spectra from the summed image 
and the sum of the individual power spectra are 
given in Fig. 0a) and Fig. 0b), respectively. In 
contrast to Fig. 0 the Thon rings are now stronger 
in the power spectrum of the sum of the frames 
with the relative strengths being essentially what is 
expected from 25 images of the same object. 

In the second control experiment, a series of 25 
images like that in Fig.[2]were taken at adjacent but 
non-overlapping areas of the pre-irradiated carbon. 
As there was a slight variation in the height of the 
carbon film at the different locations the objective 
lens current was adjusted in order to keep the vari- 
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(a) (b) 



Figure 3: Comparison from a composite exposure of (a) the 
NWPS obtained from the sum of all the frames and (b) the 
sum of the NWPS of the individual frames. The composite 
image consisted of 25 single frames taken from different dose 
fractionated images of non-overlapping areas in a region of 
pre-irradiated amorphous carbon. The imaging conditions 
were the same as in Fig.[2]with the grey scales of the images 
set as in Fig.^ 

ation in the defocus of the images to within ±90 A. 
A composite image of 25 frames was then generated 
by taking one frame from each of the 25 images. As 
in Figs. ® and [2] the corresponding power spectra 
from the Fourier transform of the sum of the frames 
from the composite image is given in Fig. [3f a) while 
that of the sum of the power spectra of the indi¬ 
vidual frames in the composite image is shown in 
Fig.§b). Like Fig. ® the Thon rings in Fig. [3] are 
more visible in the sum of the power spectra than 
in the power spectrum of the summed image. As 
expected Fig. ®(b) and Fig. ®(b) are almost identi¬ 
cal since they both consist of the sum of 25 power 
spectra from images of areas with the same thick¬ 
ness of stabilised carbon taken with essentially the 
same defocus and number of electrons. 

In Fig. ® the carbon atom positions are effec¬ 
tively the same in all the frames while in Fig. [3] 
there no correlations between atom positions in the 
different frames. The similarity of Fig. ® to Fig. [3] 
indicates that water molecules in amorphous ice are 
moving to uncorrelated positions during an expo¬ 
sure. 



Figure 4: Circularly averaged power spectra, Wi 20 ,m(u) as 
a function of spatial frequency, u, from 120 frames of image 
#230804 used in Fig. ^showing the behaviour with m. Re¬ 
sults for m = 1 (red), m = 2 (orange), m = 3 (green), m = 8 
(blue), m = 120 (grey) are shown. Exposure #230804 had 
2.33 e — /pixel/frame and using Eq. § with a DQE of 0.5 for 
the Falcon II gives an expected average background T = 550. 
The expected noise in the limit of zero spatial frequency from 
Eq. (§ is 50-v/rn but this decreases with the square root of 
the spatial frequency as more values are included in the cir¬ 
cular average. The vertical dashed line indicates the position 
of 1/3.7 A data used in Fig. [ 5 ] 

Fig- El shows the circularly averaged values of 
WM,m(u) at selected values of m obtained from the 
dose fractionated image series of amorphous ice in 
exposure #230804 used in Fig. ® For simplicity in 
factorisation, only the first 120 of 141 frames were 
used. The average background value, the magni¬ 
tude of the associated noise and their variation with 
m are in agreement with the predictions of Eq. (JsJ) 
and Eq. There are some residual correlations 
within, and between, frames that lead to small sys¬ 
tematic variations in WM,m(u) with to. These were 
removed by adding a small to dependent shift for 
each to > 1 to ensure that at the Nyquist frequency 
(in this case 1/2.08A) WM,m(u) = Wm,i(u )• The 
maximum shift required was for to = 120 and cor¬ 
responded to 5% of the background. 

Fig.® shows the circularly averaged values of 
WM,m{u) as a function of to at u = 1/3.7A (indi¬ 
cated by the vertical dashed line in Fig. [4]) from 
120 frames of exposure #23804. Also shown in 
Fig.® are results from a dose fractionated image 
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Figure 5: Measured Thon ring modulation from the cir¬ 
cularly averaged Wm m( u ) at u = 1/3.7A as a function 
of the number of electrons per unit area in the summed 
frames. The results shown using black bullet (•) are 
from the first 120 frames of exposure #230804 with 2.33 
e — /pixel/frame. The results shown using red bullet (•) are 
from the first 300 frames of exposure #171528 taken with the 
same magnification and defoucs but over 15 seconds using 
0.85 e — /pixel/frame on a different amporphous ice sample. 
The red and black dotted lines show fits to the correspond¬ 
ing experimental data using Eq. 0 . The solid blue curve 
lies almost on top of the dotted red curve and is simply 0.87 
times the dotted black curve. The error bars are based on 
circularly averaged estimates using Eq. §• 


(#171524) consisting of 300 frames from of a dif¬ 
ferent amorphous ice sample using the same mag¬ 
nification as image #23804 but with an exposure 
rate of 0.85 e“/pixel/frame. The dotted lines in 
Fig.[5]are fits to the measured results using Eq. |l7| ) 
and two adjustable parameters: er 2 and a scale fac¬ 
tor. The scale factor includes the electron optical 
terms and is proportional to the total number of 
electrons in the exposure, the scattering strength 
and thickness of the amorphous ice sample. Fits to 
successive dose fractioned images of a fixed area of 
amorphous ice give essentially the same value for 
CTq but the scale factor decreases as the sample is 
thinned by radiolysis. The fitted values of <Tq shown 
in Fig. [5] for images #230804 and #171524 are 0.38 
and 0.35 A 2 /(e“/A 2 ), respectively. There are many 
sources of systematic error in this analysis and at 
best these values for <Tq should be considered as 
estimates. Based on these estimates the radiation 


damage resulting from one incident 300 keV e _ /A 2 
is expected to induce a total mean squared motion 
of ~ 1.1 A 2 in the water molecules of the sample. 
For a dose fractionated exposure of M frames 


the Thon ring signal given by Eq. (18) (and illus¬ 
trated in Fig. [5j initially increases linearly with to 
but plateaus at higher m. As the corresponding 
noise, given by Eq. ([ 9 ]), grows as yfm there is an 
optimal to, and hence dose per image, with which 


to observe Thon rings. Dividing Eq. (18) by yfz and 


setting the derivative with respect to z to zero gives 
a transcendental equation with zero at 2.149. For a 
given resolution, u, the optimal number of electrons 
per unit area, d op t in a frame for observing Thon 
rings is 

d opt = 2.149/27r 2 cr^ 2 . (19) 

The peak is not very sharp but for a given dose the 
Thon rings around 1/3.7A. from amorphous ice (as¬ 
suming Oq = 0.37A 2 /(e /A 2 )) will have the high¬ 
est signal to noise ratio if each frame has, or succes¬ 
sive frames are grouped into blocks with an average 
exposure of ~ 4.0e _ /A 2 . 

Finally, while Fig. [5] gives the amplitude of the 
observed Thon ring modulation in Wm, m(u) as m. 
(and hence the exposure per image) is varied, it is 
also useful to look at how the corresponding visi¬ 
bility of the Thon ring modulation, i.e., the signal- 
to-noise ratio, varies. Fig. [6] shows a log-log plot 
of the variation in the signal-to-noise relative to 
the corresponding case with to = 1 in the noise 
whitened power spectra as a function of to. Results 
are shown for both the amorphous ice exposures 
shown in Fig. [ 5 ] as well as the pre-irradiation car¬ 
bon cases shown in Fig. [2] and [3] The results for 
amorphous ice fall between two extremes, given by 
a line varying as yfm for the case where the im¬ 
ages are of the same object in each frame (such as 
in Fig. [2] for the images of the same area of pre¬ 
irradiated carbon), and in the other extreme of a 
line varying as 1/ yfm for the case where the images 
are of uncorrelated objects in the different frames 
(such as in Fig. [ 3 ] made up of images of different 
areas of amorphous carbon). The #ro and 1/i/ro 
limits indicated by the dotted lines in Fi g. [6] simply 
follow from noting Wm ,m as defined in (171, varies 
as to in the limit where there is no movement be¬ 
tween frames but goes to a constant when there is 
no correlation between frames while the noise given 
by Eq. varies as yfm. The signal-to-noise ra¬ 
tio curves for both exposure #23804 and #151524 
initially vary as yfm but make the transition to 











Number of frames in average ( m ) 


Figure 6: Relative signal-to-noise ratio of Thon rings in the 
sum of noise whitened power spectra, Wm tti/Wm, i, from 
a dose fractionated exposure of M frames as a function of 
the number of successive images in each power spectrum, m. 
Results at u = 1/9.0 A from 24 frames of the pre-irradiated 
carbon images used for Fig.|2](T) and Fig. [ 3 ] (A) are shown. 
Results at u = 1/3.7 A from 120 frames of image #23804 
(•) with 2.33 e — /pixel and 300 frames of image #171524 (•) 
with 0.85 e — /pixel of Fig. [ 5 ] are shown. The results using 
Eq. |9j for the noise and Eq. ( |18| l) for the signal with Oq = 
0.37A 5 /e — are shown as solid lines. 


l/yfm for large m. As expected from Eq. (19) in 
going from behaving as 1 Jm to 1 /y/m the amor¬ 
phous ice curves reach a peak at m corresponding 
to ~ 4e“/A 2 . The solid lines in Fig.[6]show the pre¬ 
dicted behaviour of g(z)/yfz in which g(z) is defined 
in Eq. ( fl8] ) with z = 2'K 2 crQU 2 md at u = 1/3.7A 
with Oq = 0.37A 2 /(e _ /A 2 ) and d = 2.33 e _ /frame, 
(black —), or d = 0.85 e“/frame (red —). 


5. Discussion 

The recent introduction of higher performance 
electron detectors has allowed structures of biologi¬ 
cal macromolecules to be obtained by single particle 
cryoEM to higher resolution, with fewer particles 
and more easily than before @J. However, these 
reconstructions still require averaging substantially 
more particle images than predicted to be neces¬ 
sary by theory |15j . Charging of the specimen and 
its surroundings as well as physical movement of the 
specimen during exposure to the first few e _ /A 2 are 


two reasons for this discrepancy. 

To observe well defined Thon rings from images 
of amorphous ice the sample must be neither too 
thick nor too thin. If the ice film is too thin, the 
amplitude of the Thon ring signal will simply be too 
small to be seen amongst the background noise. If 
the ice film is too thick the Thon signal will also 
become weaker, especially at high resolution. Part 
of this reduction comes from the variation in de¬ 
focus as images of molecules at different heights 
have Thon ring zeroes at different spatial resolu¬ 
tions, and part is due to increased multiple elastic 
and inelastic scattering of the electron beam. If an 
ice film is too thick it is possible to use radiolysis to 
thin the film to the desired thickness by recording 
a succession of images. With 300 keV electrons the 
thickness of an ice film is typically reduced by 100 A 
for every 170e _ /A 2 of electron dose [16]. From the 
number of successive images required to completely 
burn through the ice film used in Fig. |T]we estimate 
that the image used in Fig. [I] was taken when the 
ice thickness was between 1000 and 1500 A. 

The high number of electrons used in Fig. [T] 
(304 e _ /A 2 ) ensures the Thon rings are clearly vis¬ 
ible but is not necessary. With the correct condi¬ 
tions Thon rings can easily seen with 26 e _ /A 2 . 
In fact the Thon ring visibility in the sum of the 
power spectra from 12 frames of amorphous ice ex¬ 
posure from Fig. |T] is comparable to that seen in 
the sum of the power spectra of same number of 
frames from the amorphous carbon results of Fig. [3] 
(see supplementary Sec. SVIII). The pattern of the 
Thon rings from amorphous ice differs from that of 
carbon, with the amorphous ice signal, as in Fig. [T] 
strongest around 1/3.7A while the corresponding 
pattern from amorphous carbon is strongest at low 
spatial frequencies. The 3.7 A resolution signal in 
amorphous ice comes from planes of next nearest 
neighbour oxygen atoms of tetrahedrally coordinated 
water molecules and the strength of this signal in¬ 
dicates that the ice rules of Bernal and Fowler m 
are applicable in amorphous ice. 

The major origin of beam induced motion in an 
amorphous ice film is from radiolysis rather than 
knock-on damage EE3- The rate of mass loss from 
a thin films is observed to be almost independent 
of film thickness and so most of the radiolysis frag¬ 
ments remain within the film PH- This so called 
cage-effect, allows many of the radiolysis fragments 
in ice to recombine but some such as O 2 and H 2 C> 2 , 
will accumulate and re-orient the surrounding wa¬ 
ter molecules. The successful fits shown in Fig. [5] 
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do however support the simple Gaussian Markov 
Process description of Sec 2.5 In practice irradia¬ 
tion of a plunge frozen amorphous ice film will also 
result in collective motion of local regions due to 
stress relief within the film and charging from the 
emission of secondary electrons. 

From our experimental observation we deduce 
that water molecules, each with a molecular mass 
of 18 Da, move an RMS distance of ~lA for each 
e _ / A 2 . A typical cryoEM exposure having 25 e _ / A 2 
will therefore result in a RMS displacement of water 
molecules by ~5 A. In the same way as the thermal 
motion of water molecules results in Brownian mo¬ 
tion, the beam-induced motion of water molecules 
will also be expected to displace embedded pro¬ 
tein molecules. In a 25e _ /A 2 exposure a protein 
molecule of 100 kDa, such as hexokinase, embed¬ 
ded in this film of randomly diffusing molecules 
would be expected to have an RMS displacement 
of ~1.0 A (hexokinase with a MW of 100 kDa has a 
20-40x smaller diffusion coefficient than water, and 
so should move ~4-6x less, according to the Stokes- 
Einstein equation for diffusion m in which the dif¬ 
fusion coefficient varies as 1/MW 1 / 3 ). For a protein 
molecule of 25 kDa, the beam-induced random mo¬ 
tion would be higher at ~1.5A, and for a ribosome 
of 2.5 MDa somewhat lower at ^0.7A. This addi¬ 
tional B-factor, which for hexokinase would corre¬ 
spond to ~25 A 2 , produces some extra image blur¬ 
ring on top of that resulting from other effects such 
as intrinsic radiation damage to the macromolecu- 
lar assembly. We can thus conclude that random, 
Brownian type of beam-induced motion of biologi¬ 
cal structures is unlikely to be one of the limiting 
factors in attaining high-resolution structures using 
single particle cryoEM approaches. Only for very 
small particles at resolutions of 2A or beyond is this 
type of beam-induced motion likely to be a funda¬ 
mental limitation. 

It can be argued that the observed degree of 
beam-induced motion of water molecules in amor¬ 
phous ice during an exposure is actually of posi¬ 
tive benefit in cryoEM. If water molecules moved 
much less during irradiation they would contribute 
a strong background, just as there is from carbon 
films, that would decrease the accuracy of the ori¬ 
entation determination in single particle cryoEM. If 
water molecules moved more then so too would the 
embedded macromolecules and the images of these 
would be more blurred. 
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Supplementary material: Thon rings from amorphous 
ice and implications of beam-induced Brownian motion 
in single particle electron cryo-microscopy 




Figure SI: Comparison of the spatial frequency dependence of the circularly averaged sum of power spectra and noise whitened 
power spectra from the 141 frame dose fractionated exposure (image #230804) of amorphous ice used for Fig. [l] The spatial 
frequency, w is measured as a fraction of the Nyquist frequency (1/2.08A). In (a), the sum of the power spectra of the individual 
images (5141,1 (w) and W 141 i(ttf)) along with FAf(w) using T = 684 are given. In (b), the circularly averaged power spectrum 
of the sum of the 141 frames (5i4i,i4i(tu) and Wi 4 i,i 4 i(iu)) along with VFi4i,i(2) are shown. 


SVI. Comparison of S(u) and W(u) from exposure #230804 of Figjl] 

The effect of noise whitening is illustrated in Fig.fsTja) which compares the spatial frequency dependence 
of the sum of the circularly averaged power spectrum, and noise whitened power spectrum, W± 4 ± t i, 

(using the notation of Eq. ([8])) obtained from the 141 frames used to generate Fig. 1 of the text. The 
variation with spatial frequency of Si 4 i,i(ie) is essentially that of The multiplicative factor, T, 

relating S'i 4 i i i(u;) and Af(w) is given in Eq. |[7|). The measured value for T of 684 is within 5% of that from 
Eq. (11) with a DQE(O) of ~0.5 and dose of 2.33 e _ /pixel. The noise whitening produced by division of 
S(w) by removes this variation, producing a fiat background above which the amorphous ice Thon 

ring signal around 1/3.7A, now enhanced by 1/J\f(w) ~ 24, can clearly be seen. Noise whitening can be 
applied to both the sum of the power spectra as shown in|Slja) and the power spectrum of the sum of the 
frames given in [sT|(b) . The lower noise in the sum of the power spectra provides a more stringent test for 
the noise whitening procedure; the success of this is clearly illustrated by the featureless flat background in 
Figs. 0 i of the main text. 

Residual correlations in the detector, in particular from the applied linear pixel gain correction, lead to 
a difference in the value for the sum of the power spectra versus the power spectrum of the image sum. This 
is illustrated in Fig.JsTJb) where the value of 1 ^ 141,141 (iu) is 5% greater than that of Wi 4 i,i. 

Fig. |SII| ( a) compares the spatial frequency dependence of the circularly averaged noise whitened sum of 
power spectra and that of noise whitened power spectrum of the sum of the frames. For convenience, an 
offset equal to Nyquist frequency value has been subtracted. The Thon ring signal is bigger in the power 
spectrum of the sum of the frames but the background noise is also bigger (by a factor of #141 ~ 12). In 
Fig. 1 of the main text, the Thon ring signal from the amorphous ice appears stronger in the sum of the 
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Figure SII: Comparison of the circularly averaged Thon ring signal from the dose fractionated exposure (image #230804) of 
amorphous ice as seen in the sum of the noise whitened power spectra (black) with that in the power spectrum of the sum of 
the image frames (green). The background offset, seen in Fig. |SI[ b), has been removed by subtracting the average value near 
the Nyquist frequency. The results from Fig. |Sl(b) shifted in this way are shown in (a) while in (b) the visibility of the Thon 
ring signal is shown by rescaling the power spectrum of the sum of the images in (a) so that it has the same noise level as sum 
of the power spectra curve. 


power spectra (Fig. 1(b)) rather than in the power spectrum of the sum of the images (Fig. 1(a)). This is 
because the images in Fig. 1 have been scaled so that background noise has the same amplitude. This effect 
of this is illustrated in Fig. |SII[ b) and is obtained by dividing the circularly averaged power spectrum of the 
sum frames by \/l41 so that the noise level in the two plots (green and black) is the same. 
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Figure Sill: Histogram of the 141 frames noise whitened power spectra from image #230804 of amorphous ice used in Fig. 1. 
The exponential distribution fit, —, to the measured histogram, —, is with T = 4.91. 


SVII. Statistics of noise whitened power spectra 

As discussed in the text, the distribution of the of values in a noise whitened power spectra of the 
Falcon II is completely described by its mean value F ~ n/DQE(0) in which n is the mean number of 
incident electrons per frame, through exponential distribution of Eq. (8). To illustrate this, the histogram of 
values in noise whitened power spectra of the 141 frames used in Fig. 1 is given in Fig. |SIII| The measured 
histogram is an excellent fit to an exponential distribution with parameter T = 4.91 (which differs by 5% 
from the estimate based on DQE(O) = 0.5 and n = 2.33e _ /pixel/frame). 
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Figure SIV: (a) Sum of noise whitened power spectra from the first 12 frames of image #230804 from Fig. 1. (b) Sum of noise 
whitened power spectra from the first 12 frames of the pre-irradiated carbon image used in Fig. 2. The same magnification and 
total dose (25.9 ~ 12 x 2.33/1.04 2 e“/A 2 ) were used in both (a) and (b). The Thon rings from amorphous ice around 1/3.7 A 
are stronger than those from the pre-irradiated carbon sample. 


SVIII. Comparison at 26 e /A 2 of Thon rings from amorphous ice with those from pre- 
irradiated amorphous carbon 


Fig. SIV[a) shows the Thon rings obtained from the sum of the first 12 noise-whitened power spectra 


corresponding to 25.9 e /A 2 from the 141 used in Fig. 1. For comparison Fig. SfWb) shows the correspond¬ 


ing sum of noise whitened power spectra from images of pre-irradiated carbon in Fig. 2. While the Thon 
rings from carbon are stronger at low spatial frequency, around 1/3.7A those from the amorphous ice sample 
are stronger. This is more clearly illustrated in the circular averages shown in Fig. |SV| 


Fig.JSVfa) shows the circular averages of the sum of the noise whitened power spectra shown in Fig. SIV 


In Fig. |SVp>) the noise whitened power spectrum of the sum of the 12 pre-irradiated carbon frames used 
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Figure SV: Comparison of the circular averages of the noise whitened power spectra obtained with ~ 26e — /A 2 from amorphous 
ice and pre-irradiated carbon. The W 12 1 (u) circular averages of the amorphous ice (blue) and pre-irradiated carbon (red) sums 
of noise whitened power spectra shown in (a). In (b) the corresponding circular average of the power spectrum, W) 2.12 
from the sum of the 12 frames of pre-irradiated amorphous carbon is shown. 


in Fig. SV is compared with the sum of 12 power spectra given in Fig. |SV[ a). As the pre-irradiated carbon 
does not move during an exposure the Thon ring signal adds coherently with increasing number of frames. 
The Thon ring modulation in the sum of the images, Wi 2 t i 2 (u), is 12 times that in the corresponding sum 
of power spectra, W\ 2 ,\{u). At the Nyquist frequency the image consists almost entirely of noise and the 
power spectra of the sum of the images is equal to the sum of the power spectra. The noise in the power 
spectrum of the sum is however \Zl2 times that in the sum of the power spectra. 
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SIX. Calculation of normalised noise power spectrum, A/"(u) 

The energy lost by incident high energy electrons passing through the lightly doped semi-conducting 
sensitive layer of a CMOS sensor such as the Falcon II, produces electron/hole pair excitations. The sensitive 
layer is surround by heavily P-doped layers into which the holes are free to escape but the negatively charged 
electron excitations are trapped due to the potential resulting from the differential doping. An incident high 
energy electron is recorded via the voltage drop on reverse biased diodes formed by a heavily N-doped 
implants on the surface of the sensitive layer resulting for collection of the electron excitations. The number 
of electron-hole pairs generated by an incident high energy electron, and hence the signal, depends on the 
length of the incident electrons trajectory through the sensitive layer. The thick sensitive layer of the 
Falcon II results in a large average signal but the thickness of the sensitive layer allows electrons to diffuse 
relatively large distances before they are finally collected. This diffusion results in a poor modulation transfer 
function, MTF, even with the relatively large (14/mi) pixel size of the Falcon II. The diffusion acts as a 
low-pass filter and while it lowers the MTF it has only a small effect on the detective quantum efficiency, 
DQE. In other words the spatial distribution of the signal deposited in the detector by an incident electron, 
described by the point spread function (PSF), is mainly determined by electron/hole diffusion and not 
scattering of the primary incident electron (at least for incident 300 keV electrons). The Fourier transform 
of the point spread function determines the intrinsic spatial frequency response of the detector and its square 
describes the noise power spectrum. In a real pixellated detector the noise power spectrum also includes 
pixel modulation terms of the form sin(a;)/a; from integration over the pixel and contributions aliased back 
from spatial frequencies beyond the Nyquist frequency of the pixellated sampling. 

As shown in )7] the A f(u) for the Falcon II detector is to a good approximation given by 

A f{u 1 v) = { sin(7r wuq) sin(7rauG) MTF 0 (\J Uq + Vq)\ (SXX) 

G ' ' 


where a is the pixel pitch, the sum is over the low order reciprocal lattice vectors needed to describe 
the aliasing, uq and vq are the spatial frequency components i.e., uq = u + G u , and MTF 0 is the Fourier 
transform of the intrinsic PSF. In [?; A f(u, v) was calculated from the MTFo obtained from the experimentally 
measured MTF. While this works well (reducing the variation in the power spectra to a few percent), the 
higher signal to noise in the longer exposures used in the present work show small systematic variations. 
The results presented in Figs. 1-3 were therefore obtained using a fit to a radial function <j)(u), effectively 
MTFo(zi) 2 , from zero to twice the Nyquist frequency of the measured noise power spectra using the functional 
form in Eq. ( SXX I. 

The simplest way to calculate A f(u, v) is to average the noise power spectra of many frames and apply 
the symmetry from the square pixel array. This leaves some residual noise and the problem of finding the 
overall normalisation. In the present work a radial function, (j>(u), was used to fit the symmetrized sum of 
the noise power spectra via 


S(u, v ) ~ sin 2 (7rau,G) sin 2 (7rauc) 4 >{\]+ v g)- 

G 


(SXXI) 


The overall normalisation was found by using a single Gaussian function in <p(u) to fit the low spatial 
frequency limit and a smooth general fit to the remainder via a b-spline fit that was penalty weighted to 
ensure smoothness and non-negative values. A more detailed description of this and its implementation will 
be presented elsewhere. 
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